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LONG-TERM  GOALS 

Terrestrial  runoff  and  river  input  dominates  urban  pollutant  loading  rates  Often  draining  direetly 
into  the  surfzone,  this  pollution  degrades  surfzone  water  quality,  leading  to  beaeh  elosures 
{e.g.,  Boehm  et  al,  2002),  increases  health  risks  {e.g.,  diarrhea  and  upper  respiratory  illness) 
{Haile  et  al,  1999),  and  eontains  both  human  viruses  {Jiang  and  Chu,  2004)  and  elevated  levels 
of  feeal  indieator  baeteria  {Reeves  et  al,  2004).  Surfzone  mixing  proeesses  disperse  and  dilute 
sueh  (and  other  types  of)  pollution.  On  smaller  length-seales  (smaller  than  the  water  depth), 
breaking-waves  and  bed-generated  turbulenee  mix  traeer.  However,  field  surfzone  observations  of 
turbulenee  previously  have  been  extremely  searee,  and  mueh  about  surfzone  small-seale 
turbulenee  is  not  known.  On  larger  seales  (10-100  m),  horizontal  dispersion  is  driven  by  surfzone 
eddies  and  meanders  assoeiated  with  shear  waves  {Oltman-Shay  et  al,  1989)  or  finite  breaking 
erest  length  {Peregrine,  1998).  Understanding  the  small  and  large  length-seale  mixing  proeesses 
important  to  predieting  the  fate  (transport,  dispersal,  and  dilution)  of  surfzone  traeers  whether 
pollution,  bacteria,  larvae,  or  nutrients. 

OBJECTIVES 

The  seientifie  objeetive  is  to  improve  understanding  and  modeling  of  dispersion  of  traeers 
(pollution,  feeal  indieator  baeteria,  fine  sediments)  within  the  nearshore  (a  few  100  m  of  the 
shoreline)  and  espeeially  within  the  surfzone  where  breaking  waves  intensify  mixing  proeesses 
and  drive  strong  mean  eurrents.  We  have  foeused  on  two  eomponents  of  analysis  of  the  HB06 
experiment  performed  in  Fall  of  2006.  The  first  is  studying  the  small-seale  turbulenee  in  the 
surfzone  due  to  breaking  waves  and  the  seeond,  is  larger  seale  surfzone  Lagrangian  drifter 
dispersion.  Other  efforts  inelude  analysis  of  dye  traeer  mixing  and  phytoplankton  patehiness  at 
HB06. 

APPROACH 

Surfzone  Turbulent  Dissipation  Rate: 

The  vertieal  strueture  of  turbulenee  in  the  surfzone  is  of  interest.  Both  breaking  waves  and 
near-sea-bed  shear  are  possible  sourees  of  turbulenee.  Here  a  key  turbulenee  statistie,  the 
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turbulent  dissipation  rate  e  is  estimated  from  Acoustic  Doppler  Velocimeters  observations 
following 

HB06  Drifter  Dispersion 

Two  types  of  dispersion  can  be  estimated  from  Lagrangian  drifter  data:  absolute  and  relative. 
Absolute  diffusivities  characterize  tracer  dispersion  averaged  over  many  releases  in  a  coordinate 
frame  fixed  to  the  common  release  point.  Relative  diffusivities  describe  the  spreading  of  the 
average  patch  in  a  coordinate  frame  fixed  to  the  center  of  mass  of  each  individual  patch. 
Dispersion  is  generally  calculated  by  examining  the  “spreading”  of  Lagrangian  data  (dye  or 
drifters).  The  separation  rate  of  Lagrangian  particles  {e.g.,  surfzone  drifters)  is  fundamentally 
related  to  the  spreading  of  a  tracer  patch.  As  they  disperse  differently  in  the  surfzone,  drifters  and 
dye  provide  complimentary  insights  into  surfzone  dispersion  dynamics. 

WORK  COMPLETED 


•  Spydell  et  al.  (2008)  is  now  in  press. 

•  Omand  et  al.  (2008)  has  been  revised  to  Limnology  and  Oceanography:  Methods. 

•  A  manuscript  intended  for  J.  Atmospheric  and  Oceanic  Tech.  (Feddersen)  is  almost  ready 
for  submission.  This  manuscript  deals  with  the  methods  of  analyzing  Acoustic  Doppler 
Velocimeter  data  for  estimating  the  turbulent  dissipation  rate. 

•  A  manuscript  describing  the  HB06  surfzone  turbulence  results  (see  below)  is  being 
prepared  for  publication. 

•  A  manuscript  intended  for  J.  Geophysical  Res.  (Spydell,  Feddersen,  Guza)  is  also  almost 
ready  for  submission.  This  manuscript  looks  at  surfzone  drifter  dispersion  from  Huntington 
Beach. 

RESULTS 

HB06  Experiment 

Observations  were  collected  from  15  September  to  17  October  2006  (800  hours)  at  Huntington 
Beach  CA  ,a  site  with  chronic  water  quality  problems.  A  cross-shore  transect  of  co-located 
pressure  sensors  and  acoustic  Doppler  Velocimeters  was  deployed  spanning  160  m  out  to  4  m 
mean  water  depth  (Fig.  1).  The  tide  range  was  nominally  ±1  m.  The  data  was  sampled  at  8  Hz. 
The  ADVs  sampled  between  0.5- 1.0  m  above  the  bed.  The  cross-  and  alongshore  coordinate  are  x 
and  y,  respectively.  The  mean  water  depth  is  given  by  h  and  the  vertical  coordinate  is  with 
=  0  m  at  the  sea-bed.  The  distance  below  the  mean  sea  surface  is  z'  =  h  —  z.  At  each  of  the 
frames,  hourly  estimates  of  significant  wave  height  Hgig,  cross-shore  energy  flux  F  =  Ecg,  mean 
alongshore  current  Vm,  and  turbulent  dissipation  rate  e  were  estimated. 
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Figure  1:  HB06  cross-shore  instrument  transect.  The  circled  numbers  indicate  instrumented  frames 
with  acoustic  Doppler  Velocimeters  (ADV)  and  pressure  sensors.  The  horizontal  red  dashed  lines 

indicate  the  typical  tide  range. 


Surfzone  turbulence  studies 


The  turbulent  dissipation  rate  e  was  significantly  larger  (by  a  factor  of  10)  inside  the  surfzone 
(blue  curve  first  400  hrs)  versus  seaward  of  the  surfzone  (red  curve  in  Fig.  2),  indicating  the 
importance  of  wave  breaking  to  turbulence  in  the  surfzone.  The  surfzone  e  is  related  to  the 
incoming  Hsig  (lower  panel  Fig.  2)  with  larger  waves  leading  to  larger  e.  Furthermore,  there  is 
significant  tidal  modulation  of  surfzone  e.  At  lower  tides  when  the  ADV  is  closer  to  the  surface 
(and  the  source  of  breaking  wave  turbulence),  e  is  stronger.  In  addition,  at  later  times  (hours 
500-700),  frame  3  is  alternately  within  and  seaward  of  the  surfzone  as  the  tide  goes  up  and  down. 
The  frame  3  e  varies  strongly  as  wave  breaking  is  turned  on  and  off  (see  strong  oscillations  in  blue 
curve  in  Fig.  2). 

We  seek  a  non-dimensional  scaling  for  the  dissipation  rate  e  in  the  surfzone.  Following  Terray  et 
al.  (1996),  who  developed  a  scaling  for  open-ocean  whitecapping  breaking  waves,  the  dissipation 
is  non-dimensionalized  as 

(1) 

dF/dx  ^  ’ 

where  dF/dx  =  d{Ecg)/dx  is  the  cross-shore  gradient  of  the  the  incoming  wave  energy  flux 
(Ecg)  which  is  estimated  from  the  frames.  The  non-dimensional  e  has  a  consistent  relationship 
with  the  non-dimensional  distance  below  the  mean  surface  z' /Hsig  at  each  frame  (Fig.  3).  This 
relationship  can  be  modeled  as  a  power  law  relationship  as 


Hs.ge 

dF  /  dx  y  Hsig  ) 


(2) 


where  C  and  7  are  fit  at  each  frame. 


At  the  3  surfzone  frames,  the  fit  7  vary  between  -1.5  and  -2  with  fit  skill  varying  between  0.4  and 
0.6  (Fig.  3).  This  best-fit  7  is  consistent  with  results  from  open-ocean  wave  breaking  which  has 


3 


Figure  2:  HB06  time  series  of  (top)  turbulent  dissipation  rate  e  at  frames  3  (blue)  and  6,  and 
(bottom)  incident  significant  wave  height  Hgig  at  frame  7. 


7  =  —1.9  (Terray  et  ah,  1996).  This  implies  that  the  mechanisms  by  which  turbulence  diffuses 
down  in  the  surfzone  are  similar  to  those  in  the  open-ocean.  However,  one  open  question  is  why 
the  constant  C  varies  so  much  between  the  frames.  This  is  being  investigated  further. 

HB06  Surfzone  Drifter  Dijfusivity  Parameterization 

Five  days,  Sep.  17,  Oct.  2,3,14,  and  Oct.  15,  2006,  of  GPS-tracked  drifter  data  was  collected  at 
Huntington  Beach,  CA  during  the  HB06  experiment.  The  drifter  data  consists  of  along-  and 
cross-shore  positions  as  a  function  of  time,  sampled  at  1  Hz.  There  was  a  total  of  20.2,  26.7,  16, 
20.5,  and  20  hours  on  drifter  data  on  each  day.  Position  errors  are  approximately  ±2  m  and  are 
correlated  in  time.  For  a  complete  description  of  the  drifters  and  data  quality  see  Schmidt  et  al. 
(2003)  and  a  thorough  description  of  the  data  processing  and  diffusivity  estimation  appears  in 
Spydell  et  al.  (2007).  The  cross-  and  alongshore  coordinates  are  x  =  0  and  y,  respectively. 

A  preliminary  attempt  to  parameterize  the  effect  of  varying  Eulerian  surfzone  conditions  on 
surfzone  mixing  during  HB06  is  made  by  examining  the  effect  of  surf  conditions  on  the 
asymptotic  diffusivity  for  each  of  the  five  days  of  Lagrangian  data  collected  at  Huntington  Beach. 
A  parameterization  linking  the  cross-shore  diffusivity  to  the  significant  wave-height  H^ig  and 
mean  period  T  is  given  in  Inman  et  al.  (1971)  with 
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Figure  3:  Nondimensional  dissipation  rate  Hsige/ {dF / dx)  against  non-dimensional  distance  below 
the  mean  sea-surface  z'/Hsig  at  frames  3,  4,  and  5  (see  legend)  when  in  the  surfzone.  The  solid 
lines  are  the  best  fit  to  a  power  law  relationship  H^ige/ (dF/dx)  =  a{z' / Hsig)'^  where  7  —1.9  at 

all  the  frames. 


for  the  mixing  of  a  passive  tracer  due  to  broken  wave  bores  in  the  surfzone.  Recently,  this 
relationship  was  obtained  for  a  simple  model  of  surfzone  cross-shore  tracer  diffusion  by  bores 
{Feddersen,  2007).  Specifically,  Feddersen  (2007)  showed  that  Kxx  in  (3)  is  the  average 
cross-shore  tracer  diffusivity  resulting  from  bores  of  height  and  period  T.  As  this  relationship 
parameterizes  the  effect  of  broken  wave  bores  entraining  and  mixing  tracers,  there  is  no  reason  to 
expect  this  relationship  to  hold  for  Lagrangian  drifters  which  (by  design)  are  not  entrained  in 
bores.  Thus,  it  is  expected  that  a  “cloud”  of  tracer  would  grow  much  more  than  an  equally  sized 
“cloud”  of  drifters  after  the  passing  of  a  single  bore.  Moreover,  the  mixing  due  to  bores  occurs  on 
wave  period  time-scales,  much  shorter  than  the  low-frequency  vortical  motions  which  disperse 
drifters  {Spydell  and  Feddersen,  2008). 

Using  the  wave  statistics  from  the  most  offshore  sensor,  the  asymptotic  cross-shore  diffusivity 
is  well  fit  by  this  relationship  (see  Fig  4-left).  In  particular,  fitting  the  data  to  (3)  results  in 
a  =  20.9  ±  5  with  a  skill  of  0.90.  Although  not  used  in  the  fit,  the  day  at  Torrey  Pines  (TP04) 
with  small  waves  Hf.  =  0.5  m  {Spydell  et  al,  2007),  is  also  reasonably  fit  by  this  relationship 
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Figure  4:  (LEFT)  The  HB06  asymptotie  eross-shore  diffusivity  (eireles)  aH^ /T  with  best-fit 

a  =  20.9  ±  5.  Wave  height  ifg  and  mean  period  T  are  estimated  seaward  of  the  surfzone  (~  160  m 
from  the  shoreline).  The  linear  best- fit  eonstrained  to  go  through  the  origin  (dashed  line)  has  a 
skill  of  0.90.  Also  shown  (but  not  ineluded  in  the  fit)  are  TP04  (Spydell  et  al,  2007)  data  points. 
(RIGHT)  The  HB06  asymptotie  alongshore  diffusivity  (eireles)  versus  (a)  yVmXb  where  Vm  is 
the  maximum  alongshore  current  and  Xb  is  the  surfzone  width.  The  best-fit  (dashed  lines)  results  in 
7  =  0.52  ±  0.08  with  skill=  0.93.  Error-bars  (approximately  68%  confidence  limits)  are  indicated 
by  vertical  lines.  Also  shown  (but  not  included  in  the  fit)  are  TP04  {Spydell  et  al,  2007)  data 
points.  This  figure  is  from  Spydell  et  al.  (2008)  in  preparation. 

whereas  the  second  day  with  the  much  larger  =  1.35  m  is  poorly  fit  (pluses  in  Pig.  4-EEPT). 

As  previously  mentioned,  this  relationship  is  meant  to  parameterize  the  effect  of  bores  on  tracers 
while  the  effect  of  bores  on  drifters  is  considerably  different.  Thus,  this  fit  of  the  HB06  data  to  (3) 
is  surprisingly  good  and  suggests  that  Hg  and  T  are  somehow  involved  in  controlling  the  low 
frequency  vortical  cross-shore  surfzone  motions  responsible  for  dispersion  {Spydell  and 
Feddersen,  2008),  although  the  exact  mechanisms  are  unknown. 

In  the  alongshore,  shear  wave  induced  turbulence  may  effect  the  dispersing  with  the  alongshore 
diffusivity  related  to  the  magnitude  of  the  alongshore  current.  The  following  scaling  is 
investigated 

(4) 

where  Vm  is  the  maximum  of  the  alongshore  current  measured  by  the  current  meters,  Xb  is  the 
break  point  (surfzone  width)  as  determined  by  energy  flux,  and  7  is  a  nondimensional  constant 
found  via  a  least-squares  fit.  This  scaling  is  motivated  by  mixing  length  theory  which  has  a 
characteristic  velocity  scale  (maximum  alongshore  current)  multiplied  by  a  length-scale  (the 
surfzone  width). 

The  data  is  well  fit  by  (4)  with  7  =  0.52  ±  0.08  and  a  skill  of  0.93  (Pig.  4-right).  Pitting  to  a 
constant  length  scale  L^,  rather  than  using  the  scaled  variable  surfzone  width  yxb  in  (4),  also 
results  in  a  good  fit  as  the  surfzone  width  was  nearly  the  same  on  all  days  except  one.  Thus,  the 
asymptotic  alongshore  diffusivity  is  proportional  to  the  maximum  alongshore  current. 
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IMPACT/APPLICATIONS 


Potential  impacts  include  improving  surfzone  and  nearshore  mixing  parameterizations  based 
upon  bulk  factors  such  as  wave  height,  wave  period,  bathymetry,  and  currents. 

RELATED  PROJECTS 

There  are  no  active  related  projects. 
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